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Lice are permanent, obligate ectoparasites for birds and mammals. Bovicola caprae causes hypersensitivity, irritation, 


dermatitis, anemia, lower weight gain, and lower productivity in goats. This study was conducted to investigate B. caprae 
by molecular methods based on the mitochondrial genome in the West and Northwest of Iran. A total of 1017 samples 
of chewing lice collected from ten cities in five provinces were identified using diagnostic keys. After DNA extraction 
and PCR, samples were sent for sequencing. Morphological results were consistent with molecular examinations. Nu- 
cleotide sequencing of samples isolated from different cities based on mitochondrial genome showed 100% intraspecific 
similarity. The sequences of B. caprae isolated in this study appeared in a branch next to the Canadian and Chinese 
samples in the phylogenetic tree with more than 90% similarity. The results of mitochondrial gene analysis in the present 
research showed that this fragment is useful for showing intraspecific similarity and species and genus differentiation 
of B. caprae. 
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Introduction [a 


ice are obligate ectoparasites of birds and 

mammals and are classified in the order Phthi- 
raptera. Four phthirapteran suborders are recognized: 
the chewing louse suborders Amblycera, Ischnocera, 
and Rhynchophthirina, and the sucking louse subor- 
der Anoplura [1-3]. They are hemi-metabolic insects 
and spend all stages of their lives on a host [2, 4, 5]. In 
these four suborders, there are 24 families, 303 gen- 
era, and about 5000 species of lice, of which about 
4000 species infest birds and about 1000 species infest 
mammals [6]. 

B. caprae, the goat-biting louse, has a brown- 
ish-red head and thorax, with a yellow abdomen and 
brown crossbands, and a truncated anterior margin 
of the head. Their life cycle lasts about 3-4 weeks and 
they spend all their lives on one host [4, 7], leading 
to a decrease in the quality and production of wool 
in goats [8]. Lice infestation is a major problem for 
small ruminants and causes serious damage to live- 
stock through mortality and reduced productivity, 
reproduction, fertility, and skin value in the tanning 
industry [9]. Symptoms of lice infestation in goats in- 
clude hypersensitivity to the protein in the lice saliva, 
which results in irritation and alopecia as the main 
clinical signs. Secondary infections may occur due to 
excessive scratching [4]. Lice can also cause hypogly- 
cemia, hypoproteinemia, and hypoalbuminemia fol- 
lowing the chronic loss of nutrients [10]. Moreover, 
lice infestation reduces productivity and reproduction 
in goats because of anemia and miscarriage. Reduced 
productivity in goats due to decreased weight gain and 
reproductive disorders is of considerable economic 
significance [4]. 

Lice species are conventionally identified based 
on morphological characteristics. However, accurate 
identification of species and subspecies of lice based 
on morphological characteristics is difficult because 
lice have a large variety and similar morphology. To 
solve this problem, researchers have used molecular 
markers, including mitochondrial genes and nuclear 
ribosomal genes. An appropriate genetic marker is a 
basic prerequisite for success in many evolutionary 
studies [11]. mtDNA is a valuable evolutionary tool 
for diverse structural and evolutionary aspects. These 
features include easy isolation, high copy number, no 
recombination, protected sequence and structure in 
metazoans, and a wide range of mutations in different 
molecular regions. Certain characteristics of COX1 
make it a unique and suitable marker for evolutionary 
studies. Its size and structure have been preserved in 
the studied aerobic organisms, and mutation studies 
have mapped its reaction centers, facilitating the in- 
terpretation of sequence differences in gene function 
[11]. 
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As generally agreed, lice are categorized in 
four suborders. In contrast, there is no agreement 
on the phylogenetic relationships of these groups 
and their classification [12]. Although phylogeny at 
the suborder level was suggested by Lyal in 1985, it 
was flawed due to the questionable monophyly of 
Ischnocera as a suborder [6]. Since then, many stud- 
ies have addressed phylogeny at the subordinate level 
in lice. Johnson et al. (2002) analyzed the sequence of 
three genes (EF1-18S-cox1) of 21 species from four 
suborders and showed that Ischnocera is monophy- 
letic [13]. In another study, Yoshizawa et al. (2003) 
analyzed the rrnL, rrnS sequence of 18 species and 
showed that Ischnocera was paraphyletic, grouping 
the species of Trichodectidae and Anoplura togeth- 

r [14]. Genetic analysis of insect species provides 
useful information on the taxonomic relationships, 
epidemiology, disease transmission, and control. 
Despite the importance of the economic damage of 
lice infestation, little genetic evaluation has been per- 
formed in this field. Therefore, to fill this gap, the B. 
caprae DNA sequence was analyzed using a cox1 mi- 
tochondrial marker for the first time in Iran. 


Result M 


All collected specimens were of the B. caprae 
species. Table 1 shows the total number of sam- 
ples collected per city. Using primers (HCO2198, 
LCO1490), coxl was amplified in ten samples of a 
669 bp fragment, and the PCR product was observed 
on 1.5% agarose gel. After sequencing the purified 
PCR samples of coxl, the sequences were validated 
in the BLAST system of NCBI and compared with 
the reference sequences in the GenBank. Next, they 
were registered in the GenBank with assigned access 
numbers (ID: OK135715-OK135724). The sequences 
were aligned using MEGA software and the ClustalW 
method. Phylogenetic relationships were investigated 
by drawing a phylogenetic tree based on the maxi- 
mum composite likelihood method with a bootstrap 
test (1000 replications). 

The phylogenetic tree showed a high similarity in 
the sequences of the isolates in this study. All B. caprae 
samples isolated from the studied cities were located 
in one sub-branch. Samples of B. caprae cox] isolated 
from Iranian goats in this study were located next to 
the Chinese B. caprae specimen (MF927687.1). The 
alignment results of all isolates in this study showed 
that these sequences are completely similar (Figure 1). 
During the analyses, one of the sequences related to 
the city of Bahar isolate (OK135715) was compared 
with other similar sequences in GenBank and was 
examined bioinformatically (Figure 2). The phyloge- 
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Figure 1. 


Alignment results of all isolates in this study with the closest isolate in the GenBank and an outgroup sequence. The * sign 
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Discussion [UN 


In this study, lice collected from goats in ten cit- 
ies in the West and Northwest regions of Iran were 
studied morphologically and then phylogenetically. 
All collected specimens were identified as B. caprae. 
Yakhchali et al. (2006) investigated the ectoparasites 
of sheep and goats in Northwestern Iran and identi- 
fied lice and ticks as the most common ectoparasites 
of small ruminants in this region and identified the 
species of B. caprae and Linognathus stenopsis in goats 
[15]. The species they identified in this region are con- 
sistent with the present research. 

Many studies have revealed the value of the mito- 
chondrial genome for the genetic study of population 
and intraspecific and systematic phylogeny in various 
organisms, including lice [16]. The gene composition, 
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Figure 4. 

Alignment results of amino acid sequences of the isolate in 

this study (UAR89077.1 Bovicola caprae Iran) with sequences 
extracted from the GenBank. Dots (.) in the image show the 
complete amino acid similarity with the sequence of the first line 
related to the isolate of this study. The * sign in the last line of 
each row indicates the complete similarity in that row in terms of 
similar amino acids. 


gene order, and mitochondrial genome structure are 
significantly stable across bilaterian animals [17]. The 
mitochondrial genome in metazoans is typically a cir- 
cular DNA of 13-20 kb with 36-37 genes containing 
12-13 protein-encoding genes, 2 rRNA genes, and 22 
tRNA genes [16, 18]. However, the mitochondrial ge- 
nome has an unusual structure in some species of lice 
and exists as small mini-chromosomes [16]. Among 
lice species, this fragmented structure was first found 
in the body lice Pediculus humanus corporis. The 
mitochondrial genome was then identified as small 
mini-chromosomes in some other species of lice, in- 
cluding Bovicola caprae [19-21]. The coxl gene is a 
combination of highly conserved and variable regions 
that make this mitochondrial gene a molecular marker 
particularly useful for evolutionary studies [11]. The 
phylogeny of lice at the suborder level has not been 
resolved despite decades of study. Initially, all chewing 
lice from three suborders of Ischnocera, Amblycera, 
and Rhynchophthirina were collectively referred to as 
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Phylogenetic tree of the isolate in this study (UAR89077.1 Bovicola caprae Iran) and other amino acid sequences extracted 


from the GenBank 


Mallophaga, but further studies showed that Mallo- 
phaga is paraphyletic [6]. Cruickshank et al. (2001) 
analyzed the EF1 sequence in 127 species from four 
suborders and reported that Ischnocera is paraphy- 
letic [22]. Barker et al. (2003) analyzed 18SrRNA se- 
quences from 33 species and reported Ischnocera as 
monophyletic [23]. Yoshizawa et al. (2010) analyzed 
five genes (18S, histon3, wingless, rrnL, cox1) and 
revealed Ischnocera as paraphyletic [24]. Recently, 
Johnson et al. (2018) analyzed nuclear genes of 46 
species and reported Ischnocera as paraphyletic, and 
two species of Trichodectidae were grouped with the 
species of Anoplura and Rhynchophthirina [25]. An- 
other study examined the nucleotide and amino acid 
sequence of protein-encoding genes atp6, atp8, nad2, 
nad4, nad4l, nad6, cox1, cox2, cox3, cytb, rrnS, and 
rrnL from three species of Hoplopleura. The atp8 
gene had the highest nucleotide diversity while the 
least diversity was observed in the cox] gene. The lat- 
ter finding is consistent with the present study [16]. 
Al-Shahrani et al. (2017) investigated phylogenetic 
differences based on two genes cox1 and cytb in hu- 
man head lice, and divided these lice into three class- 
es [26]. Moreover, Mokhtar et al. (2019) examined 
the genetic diversity of human head lice using the 
genetic marker coxl in Malaysia and reported that 
the collected samples belonged to clades A, B, and 
D [27]. 

This study is the first to demonstrate genet- 
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ic variation in B. caprae lice collected from Western 
and Northwestern Iran using the coxl marker. The 
results of morphological studies were consistent with 
molecular results and in general, 100% intraspecific 
similarity was observed in the nucleotide sequence 
of B. caprae samples isolated from ten cities based 
on the cox1 marker. Moreover, the alignment of the 
amino acid sequences also showed 100% similarity. B. 
caprae nucleic acid sequences in this study had 97.1% 
similarity to B. caprae (MF927687.1) from China 
and 77.1% similarity to B. bovis (MH001191.1) from 
Canada. B. caprae amino acid sequences in this study 
had 97.3% similarity to B. caprae (AUV47083) from 
China and 94.1% similarity to B. bovis (AYC65832) 
from Canada. The results of this study showed that 
cox is a useful marker to show intraspecific similar- 
ity. It is noteworthy that the coxl gene sequence in 
this research contributed to the reference sequences 
available in the GenBank and also acts as a basis for a 
larger library of a goat chewing lice sequences in the 
West and Northwest of Iran. 

The gene sequence of cox! of the samples isolated 
in the present study is partial. Therefore, a compre- 
hensive study is recommended to examine the com- 
plete sequences of this gene and compare and contrast 
comprehensive bioinformatic studies. 
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Figure 6. 
The mentioned provinces in the current study 


Materials and Methods I) 


Sampling 

This descriptive cross-sectional study was conducted during 
September 2017-March 2018 in five provinces in the West and 
Northwest of Iran (Figure 6). Samples were collected by direct 
sampling from goats’ bodies in cities Urmia, Mahabad, Sanandaj, 
Kamyaran, Kermanshah, Gilangharb, Khorramabad, Kuhdasht, 
Bahar, and Hamedan. Out of 420 examined goats, 120 animals 
were infested by lice. Collected lice samples were fixed in 70% 
ethanol and transferred to the parasitology laboratory of Urmia 
Faculty of Veterinary Medicine and were identified using valid 
identification keys [7, 28, 29]. 


Morphological identification 


All collected specimens of B. caprae were identified using 
morphological features as follows: 

e Round head and antenna with three segments 

e Distinct ocular points behind the antennae 

e The ventral surface of the thorax with dark-colored 
plates 

e All legs of similar size and with a single nail 

e Cube abdomen with black side stripes 

e Antennae in males are slightly larger than in females 
and have transverse bands 

e Each abdominal segment has a middle spot and a row of 
short hair between the groove and the spots of each segment 


DNA extraction and molecular diagnosis 


Genomic DNA was extracted using a commercial DNA 
extraction kit of MBST (MBST, Tehran, Iran) according to the 
manufacturer's instructions. The DNA quality and concentration 
of each sample were evaluated using NanoDrop (Thermo Scien- 
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tific 2000c, USA) by spectrophotometry and stored at -20°C until 
further evaluations. Primers designed by Folmer et al. (1994) as 
LCO1490 (5 - GGT CAA CAA ATC ATA AAG ATA TTG G-3) 
and HC02198 (5 -TAA ACT TCA GGG TGA CCA AAA AAT 
CA-3) were used to amplify a 669 bp fragment of cox] [30]. PCR 
was performed in a volume of 50 ul and each microtube contained 
5 ul DNA template, 5 ul 10X PCR buffer, 1 ul dNTPs (200 uM), 
4 ul MgCl2 (50 mM), 1 ul of each primer (20 uM), and 1 ul Taq 
polymerase (Sinaclon, Iran). The PCR steps consisted of an initial 
DNA denaturation stage at 95°C for 5 min, and then 35 repeti- 
tions, each cycle involving denaturation at 95°C for 45 sec, primer 
annealing at 55°C for 45 sec, extension at 72°C for 45 sec, and 
a final extension step to complete polymerization at 72°C for 10 
min. The PCR product was visualized using 1.5% agarose gel and 
UV-Transilluminator (BTS-20M, Japan). Finally, the PCR product 
was purified and sent along with forward and reverse primers to 
Takapouzist Co. (Tehran, Iran) for sequencing. 


Sequencing and genomic analysis 


Sequences were entered on the NCBI website to search for 
reference sequences with the highest similarity, and the BLAST 
method was used to find the positions of coxl. Afterwards, the 
sequences were aligned using MEGA software and any alignment 
error was resolved by the ClustalW method. All nucleotide se- 
quences obtained in the GenBank were recorded with assigned 
access numbers, and phylogenetic relationships were investigated. 
To this end, a phylogenetic tree was drawn based on the maximum 
composite likelihood method with bootstrap test (1000 replicates) 
analysis 
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